Nematodynamics with odd and rotational viscosities

We explore a novel mechanism of interactions between nematic order and flow including odd and rotational viscosities, and investigate activity-induced instabilities in the framework of this model. We show how these modes of viscous dissipation can be incorporated in the Ericksen–Leslie formalism, but it does not eliminate deficiencies of the approach based on Onsager’s reciprocal relations that lead to spurious instabilities. The suggested way of deriving nematodynamic equations, based on a specific mechanism applicable to rigid rods, is not universal, but it avoids referring to Onsager’s relations and avoids spurious instabilities in the absence of an active inputs. The model is further applied to the analysis of instabilities in active media


Introduction
The renewed interest in dynamics of nematic fluids has arisen owing to its wide applications in studies of active matter.Nematic order is commonly encountered in biological tissues [1][2][3], cells [4,5], and bacterial swarms [6,7] colonies [8], and biofilms [9].In collectives of rodlike particles [10][11][12] and biopolymers [13][14][15], nematic orientation is combined with polar activity.In all these systems, the constitutive elementary units are macroscopic, in contrast to molecular units of common nematic liquids.Nevertheless, theory of active nematic is commonly based on the established nematodynamic theory.The original Ericksen-Leslie approach [16][17][18] employing the nematic director as the order parameter, has been later extended to allow for a variable modulus through the use of the tensor order parameter [26,27].Although nematic active matter evolves far from equilibrium, its theory is commonly based on theory of passive nematics supplemented by a phenomenological active input.
The established theories aspire to be universal: they do not assume any specific mechanism of interactions between rearrangements of nematic order and flow, but derive dynamics close to equilibrium by assigning general linear relationships between thermodynamic forces and fluxes and establish relations between their coefficients via Onsager's reciprocal relations.There are general objections to this approach [28] on the basis that a combination of variables even and odd under time reversal is "physically unsound", and cannot ensure evolution to thermodynamic equilibrium.Indeed, it has been recently shown by straightforward computation a e-mail: pismen@technion.ac.il (corresponding author) [29] that instability may arise in the Ericksen-Leslie equations even in the absence of an active input, when it is ostensibly forbidden by Onsager's relations.
This paper suggests a simple specific mechanism of interaction between nematic order and flow, which is most likely to be applicable to systems with macroscopic rodlike elementary units.This mechanism generates viscous anisotropy involving both rotational and odd viscosity.The latter has been introduced in the context of superfluidity and quantum Hall effect [19] and otherwise is known to be present either in chiral media [20,21] or in fluids containing spinning components, see [22,23], as well as the review [24] and references therein.We explore how these effects modify the Ericksen-Leslie approach, but, of course, they cannot eliminate spurious instabilities.We continue with extending the model to accommodate a variable modulus.

Friction effects in a passive nematic 2.1 Antisymmetric part of the viscous tensor
Rather than describing the hydrodynamic effect of nematic rearrangement in the standard roundabout way based on Onsager's reciprocity relations, we derive nematodynamic equations directly by considering the mechanisms of both flow effect on local orientation and frictional dissipation (additional to that due to mean flow) caused by translation and rotation of oriented particles relative to the prevailing average velocity.
A velocity gradient along a rigid (rodlike) oriented particle (n • ∇)v causes the director n to rotate with the angular velocity (1) Here and further on, summation over repeating indices is implied, and 2D expressions are reduced from 3D by dropping the index a.Rotation is counteracted by the effective friction force per unit volume f a linear in ω a , which can be transformed in 3D to the stress tensor σ ij = ε ija f a .The explicit 3D form (excluding so far viscous coefficients) is where s kl = ∂ k v l is the strain tensor, which is commonly separated into symmetric and antisymmetric (vorticity) parts, s ± ij = 1 2 (s ij ± s ji ).The resulting viscous tensor η ijkl in the stress-strain relation σ ij = η ijkl s kl has a general form (excluding so far relevant coefficients) This fourth-order tensor is antisymmetric in i, j but has no definite symmetry in k, l.It can be split into symmetric antisymmetric parts, η ± ijkl = 1 2 (η ijkl ± η ijlk ), which can be assigned distinct viscosities η ± : According to the general classification of viscous tensors (Table 1 in Ref. [25]), the coefficient η − ≡ η r in the tensor η − ijkl odd in both pairs of indices i, j and k, l is identified as the rotational viscosity, and the coefficient η + ≡ η o in the tensor η + ijkl odd in i, j but even in k, l, as odd viscosity.These anisotropic terms supplement the common isotropic viscosity tensor η s ijkl = η s δ ik δ jl .

Linear stability analysis
Due to director's rotation, flow is also affected by the change of elastic energy per unit volume, defined by an appropriate Landau-de Gennes Lagrangian L. The molecular field h resulting from virtual orientational distortions is defined in the standard way [18] as which generates the distortion stress tensor Combining all these contributions brings the Navier-Stokes equation to the general form where the stress tensor σ a = σ o + σ r antisymmetric in i, j combines the odd and rotational stresses; For the purpose of linear stability analysis of Eq. ( 6) to O( ) perturbations of the ordered state in 2D oriented, say, along the axis x 1 , it is advantageous to express the perturbation velocity v through the stream function, v i = ε ij ∂ j Ψ , and expand the latter as where k is a perturbation wave vector.With the base state n 0 1 = 1, n 0 2 = 0, the indices b, l in Eq. ( 3) should be set to 1. Hence, the only O(1) elements of this tensor are , where ϕ is the angle between the perturbation direction and the unperturbed orientation.
For an incompressible fluid, taking the curl of Eq. ( 6) eliminates pressure.This brings the Fourier transform of Eq. ( 6) to The distortion stress vanishes upon linearisation, so that the influence of nematic alignment on hydrodynamics is restricted to O( 2), and linear stability of Eq. ( 6) is evident.The dynamic equation of n obtained by varying the Lagrangian is stable as well.

Modified Ericksen-Leslie approach
We explore now how adding odd and rotational viscosities would modify the traditional approach of the Ericksen-Leslie nematodynamic theory [16][17][18], which couples changes of nematic orientation and fluid motion by tracing the change of entropy in the course of evolution and requiring the resulting equations to satisfy Onsager's reciprocity relations.This subsection is purely illustrative, since odd viscosity implies broken time-reversal symmetry of the system [19], and hence, we do not expect Onsager's relations to be applicable.However, the standard approach also involves a term odd under time reversal, N = ∂ t n − s − • n, which, nevertheless, did not disqualify this approach, still considered to be standard.The standard theory starts from the equation of entropy change With a commonly used symmetric viscous tensor, (8) depends on the symmetric strain only, but odd and rotational viscosity tensors antisymmetric in i, j necessitate adding here the antisymmetric strain as well, so that the total change of entropy is presented as According to the standard procedure, the contributions to the entropy source are presented as products of fluxes by conjugate forces with coefficients having the dimension of viscosity.In the standard formalism, h is identified as the force conjugate to N and σ s ij as the force conjugate to s + ij .Including antisymmetric viscous stresses brings about additional pair: σ a ij as the force conjugate to s − ij .The forces are tied to the fluxes by linear relations with coefficients having the dimension of viscosity arranged into tensors of an appropriate rank, which, by the Onsager principle, should be symmetric to interchange of forces and fluxes.The amended relations take the general form Since the structure of the matrices in Eq. ( 12) must be compatible with the local symmetry of the nematic, they must contain only the director n and be symmetric to its reversal.Thus, the matrices S ± and B ± must contain an even number of n's.Since h and N are odd while s and σ s are even under reversal of n, the matrix G ± should contain at least a single n.The matrix S + ijkl is symmetric to the transposition of the two pairs of indices i, j and k, l, while S − ijkl is antisymmetric in i, j but has no symmetry to the transposition of k, l, since odd and rotational viscosities are not separated here.The matrix B ijkl is symmetric to the transposition of the first and antisymmetric to the transposition of the second pair of indices.The matrix G + ijk is symmetric and G − ijk is antisymmetric to the transposition of the first two indices.
To avoid proliferation of parameters, we explore a simplified version with anisotropy retained only in the matrix S − to comply with Eq. (3).With B ijkl = βδ ij ε kl , the contribution of the isotropic part of σ s to σ a vanishes, but is restored by adding the anisotropic part β n i n j δ ij ε kl .Other tensors are simplified as The dynamic equation of the director follows from Eq. ( 12): where the mobility coefficient is Γ = γ −1 , and the alignment parameters χ ± = γ ± /γ determine the response of the director to local shear and rotation.Collecting all components of the total stress leads to the enhanced nematic flow equation where σ p is the part of the total stress related to nematic alignment, usually called passive stress in the theory of active nematics, which includes both the distortion stress and the interaction terms from Eq. ( 12): The time derivative in Eqs. ( 13), ( 14) can be replaced by the substantial derivative D as in Eq. ( 6).Instabilities, originating in the terms in Eq. ( 15) additional to σ d introduced trough Onsager's relations have been shown to occur in the absence of asymmetric viscous stress [29,30].There is no chance that the described enhancements, which introduce additional parameters but lack strong stabilising factors, would be able to counteract the instability under all circumstances.This, in effect, invalidates the entire standard theory, as spontaneous deviations from equilibrium are impossible in the absence of either internal or external energy sources.

Instabilities in an active system 3.1 Instabilities in the director-based description
Evidently, instabilities may only arise in active systems.Activity is commonly introduced as an active stress σ z ij = −ζn i n j with the parameter ζ positive for extensile and negative for contractile activity.This adds to the Navier-Stokes equation ( 6) an orientationdependent term, and coupling between nematic alignment and flow may be now destabilising.As an illustration, we consider linear stability of an ordered state.The unperturbed alignment direction is irrelevant, and can be chosen along the x 1 coordinate.Due to the rotational symmetry, 2D formulation is sufficient, so that the base state is n 0 1 = 1, n 0 2 = 0. Upon linearisation and applying the operator ε il ∂ j ∂ l , the active stress supplements the Fourier-transformed hydrodynamic equation ( 7) by the term The evolution equation of the nematic director derived by varying the energy functional is supplemented by the local rotational term ε ij ωn j with ω given by Eq. ( 1).It also retains the global rotational term s − ij n j , which takes into account that the energy does not change when both the nematic alignment and the fluid as a whole are rotated with the same angular velocity.Thus, advection-alignment coupling is retained without relying on Onsager's principle, with the antisymmetric viscosity stress being responsible for rotational anisotropy.This input, specific to the interaction mechanism described in Sect.2.1, governs the flow alignment of the director, replacing the alignment parameter that appears in the Ericksen-Leslie theory, but it is not accompanied by reverse coupling generating spurious instabilities [29].Another substantial difference is the absence of a symmetric convective term s + ij n j .The 2D equation reads 50 Page 4 of 6 Eur. Phys.J. (2024) 47 :50 The molecular field h is defined in the standard way by Eq. ( 5) derived by varying the energy functional F = Ldx.The Landau-de Gennes Lagrangian is expressed in 2D as where K 1 , K 2 are splay and bend elasticities (Frank constants).This leads to the molecular field with the components We expand the perturbation in the Fourier series n i = n i (k)e ik•x dk.Upon linearisation, the indices b, j, l in the first term on the right-hand side of Eq. ( 17) are set to 1, so that i = k = 2 and only n 2 is coupled to hydrodynamics, with the product of the Levy-Civita symbols coming to −1 and ∂ 1 v 2 = −∂ 2  1 Ψ .Accordingly, the Fourier transform of Eq. ( 17) reads Stability is determined by the eigenvalues of the matrix J defining the linearised dynamics of Ψ, n 2 .Eq. ( 20) falls out, but in the case an instability develops, n 1 sets off as well, thereby changing the modulus and invalidating the director-based description.The elements of J are The trace of J is always negative, so that an oscillatory instability is excluded, and its determinant Δ is Activity dominates at long scales (k 1), far exceeding the healing length, so that the absolute value of the instability threshold |ζ * | = O(L −2 ) decreases with increasing system's size L, and vanishes in an infinite system.Since the trigonometric function in Eq. ( 23) changes sign with ϕ, both extensile and contractile activity are destabilizing in alternative ranges of perturbation angles.In a finite system, the direction of the most dangerous perturbation also depends on the ratio of splay and bend elasticities, as well as on symmetric and antisymmetric viscosities and the shape of the enclosure.

Vector-based stability analysis
A common deficiency of director-based theories is lack of conservation of the modulus |n| = 1 in Eq. ( 13), made clear by taking its dot product with n.This feature is retained in Eq. ( 17) since, generally, n • h = 0, although, unlike Eq. ( 13), other terms do vanish there.A variable modulus is commonly accounted for by taking as the order parameter a symmetric traceless matrix Q with the elements where is the modulus and d is the number of dimensions [26].In 2D, the more advantageous form of the order parameter, replacing Q, is the vector q = {cos 2θ, sin 2θ}, introduced as the nemator in Ref. [31] and independently applied by this author [29,30].The matrix with the elements be constructed by merging q and ε • q.The presence of the doubled inclination angle θ ensures the required invariance under rotation by π.
As before, we restrict to 2D, which suffices for stability analysis.Allowing for a variable modulus in vector nematodynamics based on Onsager's reciprocal relations ensures long-scale stability in a passive system, but the full expressions are quite complicated and do not exclude a short scale instability on wavelengths approaching the healing length [29].The cause of stabilisation lies in perturbations of the modulus coming with a lower degree of the wavenumber.
Similar to Sect.3.1, the basic ordered state is q 1 = = 1, q 2 = 0, but, since the change of the modulus is allowed, there are now three coupled dynamic variables, as not only q 2 but also q 1 can be perturbed.The rigid rod model with the director-based equations ( 1)-( 3) can be retained here.Rather than translating n to q, which would involve clumsy trigonometry, these equations can be used directly for the purpose of linear stability analysis, but, since the nemator rotates twice faster than the director, both flow-dependent terms are doubled.Thus, Eq. ( 17) is modified, after replacing the indices b, j, l → 1 in the first term as before, to The Landau-de Gennes Lagrangian based on the vector q is expressed in 2D as This leads to the molecular field [29,30] The algebraic term in Eq. ( 26) is linearised as −2 q 1 , so that the Fourier transform of the perturbed Eq. ( 24) is expressed as The active stress is commonly defined as σ z ij = −ζQ ij , which supplements the Fourier-transformed hydrodynamic equation ( 7) by the same term as in Eq. ( 16) with n i replaced by q i .The Jacobi matrix of the system (28) supplemented by the hydrodynamic equation is ) where H ij are the elements of the matrix  The change of the modulus does not eliminate the domination of activity at long scales in this model, as the leading term in Eq. ( 23) differs from that in Eq. ( 31) just by the factor 2. An oscillatory instability is excluded here as well.

Conclusion
The suggested mechanism incorporating odd and rotational viscosity, unusual in applications to media with nematic symmetry, may be only one of possible mechanisms responsible for interactions between flow and nematic order.The established way to derive these interactions through Onsager's reciprocal relations, aspiring to be model-independent, turns out to be faulty as it allows for spurious instabilities.This may invalidate the results of certain derivations and simulations based on the established nematodynamic theory.
It is possible that no universal theory would be ever able to capture nematodynamic interactions in all applications and on scales ranging from molecular to macroscopic.Flow-alignment interactions and the way they are affected by activity are likely to be specific in different applications, especially biologically related, and require further deep insights.Activity, defined here in the standard manner, may be also introduced in different ways corresponding to specific mechanisms.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.

Funding
Information Open access funding provided by Technion -Israel Institute of Technology.Data Availability Statement No data are associated with the manuscript.